Analysis of complete positivity conditions for quantum qutrit channels 
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^ We present an analysis of complete positivity (CP) constraints on qutrit quantum chan- 

C^l nels that have a form of affine transformations of generalized Bloch vector. For diagonal 

(~| , (damping) channels we derive conditions analogous to the ones that in qubit case produce 

■ 

^ ■ tetrahedron structure in the channel parameter space. 

■ I. INTRODUCTION In our work, we aim at presenting similar 

> , 

CN ' analysis for qutrit channels. Qutrit states are 

00 

OO ' states belonging to three dimensional Hilbert 

\ The analysis of quantum channels - corn- 

space and, in analogy to qubit case, one can 

use generalized Bloch formalism to describe 
\ of the crucial points for quantum informa- V 

• • , their evolution [4|. Generalized Bloch equa- 



■ pletely positive trace preserving maps - is one 



tion theory. Quantum channels correspond to 
processes that physically may take place and 
lead to evolution of a quantum state 
For qubit case, we already know the analysis 
of quantum qubit channels Q]. This analy- 
sis gives us a connection between the physi- 
cal evolution of the system written via Bloch 
equations and quantum channel formalism [3(. 
Hence one can derive mathematical conditions 
on parameters that appear in Bloch formalism. 
These mathematical conditions are the conse- 
quence of the fact that each physical process is 
a completely positive map. 



tions, that describe for instance three level 
atoms, may present the physical context for 
the evolution of qutrit quantum state. In anal- 
ogy to qubit case, we investigate the evolution 
of generalized Bloch vector that evolves within 
a Bloch ball. As a natural choice, we inves- 
tigate qutrit channels of the general form of 
linear transformation on the qutrit Bloch vec- 
tor - studying affine transformations on qutrit 
Bloch vectors. We parameterize qutrit chan- 
nels and then derive conditions for channel pa- 
rameters in order to obtain physical transfor- 
mations: completely positive maps (CPM). In 



'Electronic address: |agata.checinska@fuw.ed"irpll qubit channel analysis a tetrahedron structure 



of completely positive maps appears (for qubit where *-product is defined as {A * B)i = 

unital channels), we show analogous analysis dij^AjB^, with dijk being totally symmetric 

for qutrit channels, for which more sophisti- tensor [4]. Qutrit states belong to a general- 

cated channel geometry emerge. This new re- ized Bloch ball, qutrit pure states belong to 

suit on qutrit channels can be linked with the the unit sphere 5^ = {n S : = 1}. 

analysis of bipartite qutrit states via Jami- However, physical qutrit states do not occupy 

olkowski isomorphism. entirely the generalized Bloch ball. The pure 

qutrit states (states satisfying conditions ([3])) 



II. STATE DESCRIPTION form a subset of the unit sphere. 

parametrized with 4 parameters 

Let us first recall the idea behind the Bloch 



hey can be 
6(] as follows 



formalism. This, in qubit case, corresponds to _ e*^i gin 6* cos (/>|0) + e*'^^ sin6'sin(/)|l) + 

the choice of representation of the qubit state „,„, 

+ cos 6*12), (4) 

density operator: the basis of Pauli matrices cjj 

where < 6',</) < f, < xi>X2 < 27r and 
= -(1 + 6(t), (1) overall phase was omitted. Hence, the set of 



pure qutrit states is a 4 dimensional subset of 
7 dimensional sphere. 



2 

where 6 is a three dimensional, real Bloch vec- 
tor, describing the qubit state and satisfying 
iP' <1 (equality for pure states). Qubit states 

occupy entirely the Bloch bah. In a similar m. COMPLETELY POSITIVE TRACE 

way we can represent a qutrit state, a state PRESERVING MAPS 

belonging to three dimensional Hilbert space. |-| 

In qutrit case the choice of representation is set "^^^ shown m [3| that physical transfor- 

to be the basis of Gell-Mann matrices A^, the mations must not only be positivity preserving 

generators of 5*^7(3) group [6] there exist more subtle conditions to sat- 

1 _ ^ isfy, these are called complete positivity (CP) 

p,t = -{i + V^n\). (2) nh 

conditions [2, lZ[. The classification of qubit 

Here n is a generalized Bloch vector, real and , , j- j. i j. -j.- -j. ■ 

° ' channels according to complete positivity is 

eight dimensional. Qutrit states can be char- ,,1 Mi j. j. 4. ■ -i 

° ^ well known [^J. We want to present similar 

acterized by condition n < 1. However, qutrit ^ ■ c j. -j. 1 1 

— ' ^ analysis for qutrit channels. 

case is more sophisticated: pure states are 1, • i -j-u tr-iu 4- 

^ ^ We represent physical system with Hiibert 

states for which two conditions are satisfied n, ntnj\ ■ .1 11 c ^^^ 11 

1-1 1-1 space ri. o[rL) is the algebra of all bounded op- 
erators on "H, a linear map $ : B{T-i) ^ B{7i) 



= 1, n* n = n, (3) is completely positive if for every positive inte- 



3 

ger m the map: partial trace of with respect to the first 

$M = $®lM : ^ B{n)(g>M^"'\ subsystem (A) gives the miit operator for the 

second subsystem: TvaD^ = I. To evaluate 

, , ^(^\ , the entries of dynamical matrix, we need to 

IS positive (where ' is the identity operator 

,1 11 i Aim) r 1 , ■ compute the action of the channel ^ on Ejk- 

on the algebra Jvl^"'> oi mx m complex matri- 

a„, , , , . , , . , , , Once again, when its action is rewritten as 

. Clearly, this amounts to saying that 



ces 



N 



$ acts on a subsystem A of a larger Hilbert ^^.z n ro\ 

space and there is a reservoir (or subsystem B) a,T=i 

on which we act with unit operator . Here, where coefficients ^fj.u,aT characterize channel 

we do not know the dimension of the reservoir action, we see that we work with numbers, 
and therefore ^^^^ must be positive for any m. 

It was shown that every CPM has an operator- Channels versus states - Jamiolkowski 

sum or Kraus representation fl] isomorphism 

$'^^^(p) = ^/Cip/Cj, (6) On the other hand, iV^ x TV^, positive and 

i 

. , 1 ^ , . , r , , • r hermitian dynamical matrix must corre- 

with /Cj being a set of Kraus operators satisfy- 

i^ti- TT spond to a density operator acting on an A^^- 

mg Lji^^i = I- 

„ 1 i 1 ii • i r ^ dimensional Hilbert space. This correspon- 

io evaluate whether a given transformation 

, ,. ^ . , , , .,. , dence is up to the normalization factor, since 

(a linear mapj is completely positive we need 



to construct the so called dynamical (or Choi) 



TrD^ = N. Hence = jq is a proper 



p,, . ,t9 /^T- ,1 1- • density matrix that we can write as 

matrix of the size i\ x JS (IN is the dimension 

of the system of interest). We will denote the t i a / -ia ^ i -v / -i 

^ n Q P^ = j^D^ = — ^ 

dynamical matrix with 2, Lfl]. Dynamical *,i=i 

1 1 ^ 

matrix represents uniquely channel action. We _ _ _ (^{^E- ■) (9) 

denote with Ej]^ N x N matrix with 1 at po- *>i=i 

sition (j,k) and zeros elsewhere. The map ^> is The set of density operators defined by the dy- 

CPM iff namical matrices is only a subset of density ma- 

JL^ trices in A^^-dimensional Hilbert space, since 
= ^Eij)^Eij, (7) 

ij=i dynamical matrices must satisfy Tr^D^ = I. 

is positive semi-definite (Z)$ > 0). The fact that completely positive maps ^^"^^^ 

Channel $ must preserve hermiticity of den- which are represented uniquely by they dynam- 

sity matrix and therefore its dynamical matrix ical matrices, correspond to states is known as 

must be hermitian: = D^. Trace preserv- the Jamiolkowski isomorphism Q]. 

ing of the density operators means that the Therefore, when analyzing quantum qutrit 



channels we can reinterpret it as an analysis 
of two qutrit quantum states, in qutrit case 
= 3, therefore /o^ is a 9 x 9 matrix. 



IV. QUANTUM QUBIT CHANNELS 
A. Bloch equations 

To recall the qubit case analysis, we can 
start with a two-level quantum system and its 
evolution. The latter can be written by means 
of Bloch equations that are equations for com- 
ponents of Bloch vector b = {u,v,w). If we 
take, for instance, the decoherence of a two- 
level atom, these equations read 



u = ——u — Av, 
T 

i) = —jf-v + Alt + Qw, 

W = - — {W - Weq) - ^W, 



(10) 



and represent the evolution of the system. Here 
f], A are Rabi frequency and detuning respec- 
tively. ^ stand for decay rates for the atomic 
dipole {i = u, v) and decay rate of the atomic 
inversion {i = w). These equations, when put 
together, give rise to an affine transformation 
of the qubit Bloch vector that is governed by 
the parameters listed above. Any physical pro- 
cess amounts to a transformation of the qubit 
state that is already a completely positive map. 
However, not every affine transformation of the 
(qubit) Bloch vector will be a completely pos- 
itive map. 



B. Afflne transformations on qubit Bloch 
vectors 

The analysis of completely positive trace 
preserving maps on M2 (complex two dimen- 
sional matrices) has been studied extensively 

a a 

and gives the answer to the problem. 
Without loss of generality one can analyze 
qubit channels that transform qubit Bloch vec- 
tor according to 



^ib : b^b' = kH+t^^, 



(11) 



where matrix = diagjA^*, A2'', A3''} con- 
sists of damping eigenvalues A.f' and t"^^ = 
{t"^ ,tf' ,t'^) is a translation. The image of the 
set of pure states (6^ = 1, Bloch sphere) under 
such transformation is the ellipsoid 



( 



qb / + V 



Af 



Y + {- 



3 \2 



Af 



1, 

(12) 



with its center defined by t and its axes by Aj. 
The set of conditions on both A?'' and t'f can 
be found in . When we limit ourselves just 



to diagonal qubit channels (meaning f^^ = 0), 
then the set of allowed A^^ forms a tetrahedron 
structure . This structure reappears also 
in the space of two qubit states. 

V. QUANTUM QUTRIT CHANNELS 
A. Qutrit Bloch equations 

As before, we can start the analysis of trans- 
formations on qutrit quantum states with the 
analysis of a three level atom for wich we can 



write down Bloch equations. Three level atom the equations we have seen in the qubit case, 
is not the only possible physical realization In this work we use slightly different nota- 



Id, 



111 ] but it is very illustrative. The ana- tion for the qutrit vector - we already have in- 



log of the Bloch vector for the case of a three troduced qutrit Bloch vector n related to the 
level atom was in the beginning introduced as a choice of Gell-Mann matrices basis (in some 



(eight dimensional, real) coherent vector S [13], works generalized Bloch vectors are also called 
which components (denoted as u, f, ti)) were de- coherent vectors J]). These two vectors (5 and 
fined as n) are of course equivalent. 

Parameters that appear in the qutrit Bloch 



Ujk = Pjk + Pkj, 
Vjk = KPjk - Pkj), (13) 



equations have physical background, therefore 
the resulting affine transformation is a com- 



^i'^ Y /(/ + 1) pletely positive map. However, we can as 

x(pii + P22 + ••• + Pn — I Pi+i i+i), ^^11 ^sk the opposite question: given an ar- 



bitrary affine transformation on qutrit Bloch 
vector what are the conditions on its parame- 
ters which guarantee complete positivity? 



B. Affine transformations of qutrit Bloch 
vectors 



with 1 < j < A; < 3 and 1 < / < 2. Now, as an 
example of the physical system we can take a 
three level atom for which nonzero dipole mo- 
ments are between levels 1 and 2, and 2 and 3. 
The atom interacts with the electric field (two 
electromagnetic waves incident on the atom) 

and we assume that detunings are the same having in mind the question stated above, 
(Ai2 = -A23 = A). The corresponding Bloch ^-^^ ^^^^ transformations of qutrit Bloch 
equations for coherent vector S are ^^^^^^ ^^lat have a form 

ui2 = Avi2 + l3vi3, ^ - -/ A - , r 

: ni-^n =An + t, (15) 

U23 = -A7;23 - avi3, 



Ul3 = Pvi2 - aV23, 



where A = diag{Ai, Ag} consists of 8 damp- 
ing coefficients and t is an eight dimensional 
vi2 = - Ani2 - (3ui3 + 2awi , translation. The image of the set of pure states 

V23 = Au23 + auis - Pwi + V3(3w2, under this transformation is 

vi3 = -Pui2 + au23, ^ C'' = 1 (16) 

i=l 



Wi = -2aVi2 + PV23, 
W2 = 



together with the condition for *-product n 



-V^PV23, (14) ^ 

n = n 

where a, B are related to two Rabi frequencies , , ^/ + 

I — I n- — r.- n- — I 



I2I ]. These equations are a generalization of Aj Aj ' 
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On the other hand, parameters Aj, ti must sat- C. CPM conditions for diagonal channels 

isfy 

In qubit case, the action of the diagonal 
channel can be written as 

b^b' = A-?* = diag{Al\Af, Af }, 

^(A^n, + ti)^ < 1, (18) (19) 

i 

whereas for qutrits we have 
n^n=An, A = diag{ Ai, As}. (20) 



according to the requirement n''^ < 1. How- 
ever, complete positivity is a much stronger 



In both cases, we assume that the nature of Aj 
parameters is quasi-damping, hence |Aj| < 1. 
This comes from the fact that < 1 at all 
condition than condition saying that we can- ^-^^^^ therefore, the change of any initial Bloch 
not exceed value 1 for the length of Bloch vec- vector will lead to a vector within the (gener- 
tor. The latter, in qubit case, amounts only to alized) Bloch ball. Dynamical matrix Z)$ for a 
statement that the density operator must be a Qutrit channel of the form ^ must be positive 

■ 1 n J. T i- -i. u semi-definite in order to correspond to CPM. 

positive dennite operator, in qutrit case how- ^ 



ever, it is even less than that - since not every 
point within the 5^ sphere corresponds to den- 



(21) 



There are nine eigenvalues di that must be non- 
negative to satisfy positivity of D$. The first 

six eigenvalues give rise to conditions that can 

sity operator. , , 

^ be written as 

Transformation (fTSl) can be rewritten to give , a , 3/a an ^ n 

i - Ag + 2l^^4 — Asj > U, 

channel coefficients ^fj.u,aT- To construct dy- i _ — §(A4 — A5) > 

namical matrix D^, we apply the channel ac- 1 — Ag + |(A6 — A7) > 0, 

tion to i^jfc 1-^ ^(Ejis), representing it in the 1 — Ag — KAg — A7) > 0, 

basis of Gell-Mann matrices: Ejk = -j^ni^^a 1-^8 + K^i - ^2) + |(Ag - A3) > 0, 
(where a G {0, ...,8}, Aq = and n-'a can ^ 2^^'^ A2) + 2(Ag A3) _ 0. 

These conditions alone lead to the set of al- 
lowed Ai that has a polyhedron like struc- 
ture. In qubit case we have similar set of 
equations for A*^^ that define the tetrahedron 
structure. However, in qutrit case there are 



be interpreted as an analog of Bloch vector). 
We will first look at the channels that consist 
only of damping matrix (diagonal channels) 
and do not have a translation. These chan- 
nels are in fact unital, since they leave maxi- 

„ . 1 , ; 1 , /,! 11 1 three remaining inequalities (given by eigen- 

mally mixed state unchanged (they are called 



values dj, ds, dg > 0) which reveal coupling 
between all the parameters. 



bistochastic maps [2]). Later on, we will look at 
channels that include also translations of Bloch 
vector. d7(Ai, Ag) > 0, 
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dsiAi, As) > 0, The structure of the set {A<^^*^} 

d9(Ai,...,A8) > 0. (22) 

Because of their numerical complexity they are For qubit case, the allowed values of damp- 
discussed in Appendix. Matrix is hermi- ing parameters {A^^^^^}i=i,2,3 form a charac- 
tian, therefore eigenvalues ^7,8,9 must be real, teristic structure (tetrahedron, ^]). We are in- 
For some cases, three conditions ^ reduce terested in the structure that appears in qutrit 
to just two (see Appendix). All the inequali- case. The main obstacle here is the size of pa- 
ties characterize the set of ahowed {A"^^^^}, in rameter space. We have eight parameters on 
other words, channel parameters A^ for which which we impose our CPM constraints. We can 
$ is a CPM. The boundaries of the set { A'^^^^} investigate the {A^^^^} set projecting it onto 
can be computed by analyzing values of A^ sat- subspaces. It is easier to work with parameters 
isfying equations instead of inequalities given paired according to (Ai,A2), (A4,A5), (A6,A7) 
by ([21]) and ([22]). In principle, parameters A^ (this pairing refers to the form of Gell-Mann 
can be time dependant, still, conditions matrices). We can also put together (A3,A8) 
and ([22]) must be satisfied for any time t to (two diagonal Gell-Mann matrices) though not 
have a CPM. necessarily, since the equations are not sym- 
If we assume, for example, that metric in these two. 

On FigdJ-FiglUwe show projections of {A*^^^^} 

dn(t) 

— ^ = ^■ni{t), (23) onto the various subspaces in 8 dimensional 

space of parameters Ai,...,A8. The dark re- 
then time evolution of the Bloch vector n{t) is gions in these figures correspond to these com- 
given by: binations of Aj which are satisfying CP condi- 
tions. There are many such projections that 
ni{t) = e'^**nj(0). (24) can be obtained from the conditions that we 

have derived. In principle, the structure of the 

We can then identify Aj = e"^'* and conditions set of {A'-^^^} is not simply a generalization 

on Aj will impose conditions on 7j. For this of a tetrahedron. Since we have three (or two) 

type of evolution, one can write the Lindblad conditions for Aj that couple all the parame- 

equation for qutrit density operator p{t), cor- ters (in a nonlinear way), the simple polyhe- 

responding to the channel action. Some more dron type structure (emerging from inequali- 

details on relation between complete positivity ties that are linear in Aj, ([2T]) ) is altered. In 

and master equation and Lindblad operators the figures, one can see combination of almost 

one can find in jsl. [l^. \l\ . rough edges with smooth behavior elsewhere. 
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Y 




X 



FIG. 1: The dark region in the figure shows these 
values of parameters Ai = A2 = Y, — X 

that satisfy CP conditions. The qutrit channel has 
only the diagonal (damping) part. 



Y 00 




X 



FIG. 2: The dark region in the figure shows these 
values of parameters A3 = Ag = F, A^^s g = X 
that satisfy CP conditions. The qutrit channel has 
only the diagonal (damping) part. 

D. CPM conditions for channels based 
only on translations 

In this section we will analyze shortly the 
constraints of complete positivity on the pos- 
sible translations. The change of qutrit Bloch 
vector in this case will be of the form: 

fi 1-^ fi' = n + t, (25) 

where t = {ti,...,ts). This type of channel is 
nonunital. Below we show some examples of 



Y 00 




X 



FIG. 3: The dark region in the figure shows these 
values of parameters A3 = X, Ag, — Y, A^^s.g — 
XY that satisfy CP conditions. The qutrit channel 
has only the diagonal (damping) part. 



X 

FIG. 4: The dark region in the figure shows these 
values of parameters Ai = A2 = X, A3 = Ag = 
XY, Ai^i_2.3,8 = Y that satisfy CP conditions. The 
qutrit channel has only the diagonal (damping) 
part. 

channels, for which we choose just two free pa- 
rameters. First, let us look at the translation 
of the form 

Ti = (X,X,y, 0,0, 0,0,0). (26) 

It turns out that effectively, parameters X, Y 
must satisfy 

1 - SVSV^X^ + y2 > 0, (27) 

what graphically is represented on Figl5]- the 
dark region corresponding to CP-allowed pa- 



rameter values has an ellipsoid form. 



X, Y must satisfy 




1 + 3Y- 3\/3^2X2 + y2 > 0. (29) 

This is shown on Fig (61 The ellipsoid shape 
from FiglUis now reshaped. 

E. CPM conditions for diagonal channels 
with translations 



FIG. 5: The dark region shows these values of pa- 
rameters X, Y that satisfy CP conditions, when 
translation acting on qutrit Bloch vector has the 
form Ti = {X, X, F, 0, 0, 0, 0, 0). 



Y 00- 



V 



-02 00 



FIG. 6: The dark region shows these values of pa- 
rameters X, Y that satisfy CP conditions, when 
translation acting on qutrit Bloch vector has the 
form Ta = [X, X, F, 0, 0, 0, 0, Y). 

On the other hand, if we let the translation 
to shift also the 8th component by the same 
amount as the 3rd component, therefore trans- 
lation having a form 



r2 = (x,x,y,o,o,o,o,y). 



(28) 



then the allowed set of parameters X, Y is fur- 
ther limited with respect to X - this parameter 
must satisfy Y < 1/6. And both parameters 



We have seen what are the CP conditions 
for diagonal channels and investigated some ex- 
amples of channels built up only with trans- 
lations. What occurs when these two effects 
combine? Let us take the channel that changes 
the Bloch vector according to 



n 



n = An + t, 



p ^ <^{p) = i(l + ^(An + t) • A). (30) 

In this case, to evaluate positivity of D$ we 
evaluated the principal minors of the matrix. 
In principle, the matrix is positive when all 
the principal minors are positive, and it is neg- 
ative, when the principle minors have alter- 
nating signs. However, it may happen that 
for some parameter values the principal minors 
equal to zero and the method do not detect all 
possible parametrization allowed by CP condi- 
tions. Nevertheless, we use this method to an- 
alyze some cases and detect possible regions of 
positivity of dynamical matrix D^. We do not 
present here the list of inequalities correspond- 
ing to CP conditions because their complexity 
would unable any insight into the problem. We 
project the set of {Aj, tj}'-^^^^ onto some sub- 
spaces to gain a geometrical picture. Below, 
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we can see two examples. On FiglT] the dark 
region corresponds to the CP-ahowed values of 
X, Y, for which we assume that all Ai = X 
and ti = Y. Fig (8] shows a different choice of 



Y 00 



FIG. 7: The dark region in the figure shows these 
values of parameters Ai — X, ti = Y that satisfy 
CP conditions. The qutrit channel contains both 
the diagonal (damping) part and the translation. 




FIG. 8: The dark region in the figure shows these 
values of parameters A12 — X, Aj^i,2 — F, U ~ 
XY that satisfy CP conditions. The qutrit channel 
contains both the diagonal (damping) part and the 
translation. 

parametrization and relation between damping 
parameters and translation. We let Ai^2 have 
independent value (X) from the rest of Aj (Y), 
and we assume that translation is equal to the 
product of these two (XY) . 



F. Two-qutrit states and afltine 
transformations of qutrit Bloch vectors 

As already said, the dynamical matrix 
corresponds to a density matrix via = 
jjD^. The latter, in our case {N = 3) de- 
scribes a class of two-qutrit states that can be 
parameterized by {Aj, ti}'-^^^ . The two qutrit 
state space is being investigated. Especially, 
the so called magic simplex which can be con- 
sidered an analogof the magic tetrahedron of 



bipartite qubits 



15( 1 . The magic simplex of bi- 



partite qutrits is only embedded in the space of 
all bipartite qutrits. As an example, it does not 
contain a state given by the density operator 
1 



P 



|^)(*| + 2|$)($|), 



(31) 



where 1^-) = |0,0) and {(j)) = ^(|1,1> + 
1 2, 2)). Interestingly, this state can be obtained 
from p$ = -^D^ (N=3) with a proper choice 
of parameters: As^ej.s = 1 and the rest equal 
to 0. A diagonal channel with such parameter 
values will transform any qutrit Bloch vector 
according to 



n 



n 



{0,0,n3,0,0,n6,n7,n8}. (32) 



Geometrically, the channel projects the Bloch 
vector onto the 3-6-7-8 subspace and the other 
components of n are lost. It resembles there- 
fore a phase flip type channel [l|. 

VI. SUMMARY 

Basing on the generalized Bloch formalism 
for qutrit quantum states we have analyzed 



quantum qutrit channels that have a form of 
affine transformations on Bloch vectors. The 
aim was to derive complete positivity condi- 
tions on channel parameters that may appear 
in equations of evolution of the Bloch vec- 
tor. We analyzed diagonal channels (only with 
damping coefficients), for which we obtained 
CP conditions on parameters in form of in- 
equalities. Analogous inequalities appear in 
qubit case - for which CP-allowed channel pa- 
rameters form a tetrahedron structure. The 
structure of the corresponding set in qutrit 
case is more sophisticated, and reveals not 
only polyhedron like characteristics. We an- 
alyzed also channels which allow only shifting 
the Bloch vector, in which case we investigated 
some specific examples of translations and the 
corresponding CP constraints. The combined 
effect of damping and shifting qutrit Bloch vec- 
tor (diagonal channels with translation) was 
also presented by projecting the CP-allowed set 
of channel parameters onto specific subspaces. 
At the end we looked at the two qutrit states 
that correspond to qutrit channels we investi- 
gated, via Jamiolkowski isomorphism. As an 
example we give a channel that corresponds to 
a state which does not belong to the so called 
magic simplex. 

One of the interesting points would be to estab- 
lish the relation between the set of two qutrit 
states given by the dynamical matrix that we 
analyze and the magic simplex for qutrits. Also 
the analysis of the structure of CP-allowed 
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channel parameters with respect to entangle- 
ment breaking properties could reveal some in- 
triguing results. 
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APPENDIX 

In the discussion of the diagonal (damping) 
qutrit channels we presented a set of condi- 
tions (j21l22p that come from imposing posi- 
tivity condition on They correspond to 
eigenvalues of the matrix, and three of them 
have sophisticated form. In qutrit case, the 
dynamical matrix is nine dimensional. We 
already showed 6 eigenvalues. The remain- 
ing three correspond to finding roots of poly- 
nomial of the 3rd order. Since D$ is hermi- 
tian, all its eigenvalues must be real, to guar- 
antee complete positivity of the channel, they 
must also be nonnegative. The polynomial 
P{x) = Ax^ + Bx^ + Cx + D which we analyze 
in order to obtain the rest of CP conditions 
(eigenvalues ^7,8,9) has coefficients 

A = S, B = -24(1 + A3 + Ag), 

C = 18((A3 + Ag)'' - (Ai + Aa)^ - {A4 + As)^ - (As + Ay)^) + 

+24(1 + 2A3 + 2A8 + A3A8), 
D = -8 - 18{{A3 + Ag)^ - (Ai + Aa)^ - (A4 + Aa)^ + 

-(As + Ay)^) + 27A3((A4 + As)^' + (Ag + A-,f) + 
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-54(Ai + A2)(A4 + AsXAe + A7) + 
+9As{4{Ai +A2f + {A4 + As)^ + (Ag + A^f) + 
-24(A8 + A3 + A3A8) - 4A8{A3 + As)^ + 
-32A8A3 - 20A|A3. 

The roots can be of course found explicitly, 
but we do not present them here. For some 
specific parameter values, number of real roots 
can be reduced (and therefore number of CP 
inequalities). To deduce that one has to ana- 
lyze function 

f{a,b,c) = 3-^{3b-a'^f+54-^{9ab-27c-2a^f, 



and evaluate it at f{B/8,C/8,D/8) (where 
B,C,D are polynomial coefficients given 
above). For the dynamical matrix 
this function is always nonpositive (and 
indicates therefore real roots). When 
f{B/8,C/8,D/8) = the polynomial has 
three real roots and at least two are equal 
- then the number of CP conditions reduces. 
This occurs for parametrization shown in Figl2j 
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